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■ We present essentially exact solutions of the Schrodinger equation for three fermions in two dif- 

ferent spin states with zero-range s-wave interactions under harmonic confinement. Our approach 
. covers spherically symmetric, strictly two-dimensional, strictly one-dimensional, cigar-shaped, and 

pancake-shaped traps. In particular, we discuss the transition from quasi-one-dimensional to strictly 
one-dimensional and from quasi-two-dimensional to strictly two-dimensional geometries. We deter- 
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mine and interpret the eigenenergies of the system as a function of the trap geometry and the 



■ strength of the zero-range interactions. The eigenenergies are used to investigate the dependence 



of the second- and third-order virial coefficients, which play an important role in the virial expan- 
sion of the thermodynamic potential, on the geometry of the trap. We show that the second- and 
third-order virial coefficients for anisotropic confinement geometries are, for experimentally relevant 
C/j ■ temperatures, very well approximated by those for the spherically symmetric confinement for all 

s-wave scattering lengths. 



PACS numbers: 



I. INTRODUCTION 



\ There has been extensive interest in ultracold atom physics in the last decade [J-O- Ultracold atomic bosonic and 
fermionic gases are realized experimentally under varying external confinements. In these experiments, the number 
of particles and the scattering length of the two-body interactions are tunable [4||. Although the complete energy 
spectrum of the many-body system cannot, in general, be obtained from first principles, the energy spectra of selected 
few-body systems can, in some cases, be determined within a microscopic quantum mechanical framework [sl-fT^. In 
^ , some cases, the properties of the few-body system have then been used to predict the properties of the corresponding 
many-body system [ll| - [l9| . 

The behavior of atomic and molecular systems depends strongly on the dimensionality of the system (20l - [23j . In 
three dimensions, e.g., weakly-bound two-body s-wave states exist when the s-wave scattering length is large and 
\^ ' positive but not when it is negative. In strictly one- and two-dimensional geometries, in contrast, s-wave bound states 
. exist for all values of the s-wave scattering length Q . 

■ In ultracold atomic gases, the de Broglie wavelength of the atoms is much larger than the van der Waals length 
' that characterizes the two-body interactions. This allows one to replace the van der Waals interaction potential 

^) ■ in free-space low-energy scattering calculations by a zero-range s-wave pseudopotential j24H26l] . If the particles are 
I \ placed in an external trap, the validity of the pseudopotential treatment (at least if implemented without accounting 

■ for the energy-dependence of the coupling strength) requires that the van der Waals length is much smaller than 
° the characteristic trap length [i^,!!^. In many cases, the use of pseudopotentials greatly simplifies the theoretical 

J> . treatment. For example, the eigenequation for two particles interacting through a s-wave pseudopotential under 
r ■ harmonic confinement has been derived analytically for spherically symmetric, strictly one-dimensional, strictly two- 
r> , dimensional and anisotropic harmonic potentials [5-7]. 

■ The s-wave pseudopotential has also been applied successfully to a wide range of three-body problems, either in 
■ " " ' free space or under confinement (29l - l35j . The present paper develops an efficient numerical framework for treating 

the three-body system under anisotropic harmonic confinement. The developed formalism allows us to study the 
dependence of the three-body properties on the dimensionality of the system. We focus on fermionic systems con- 
sisting of two identical spin-up atoms and one spin-down atom. The dimensional crossover of two-component Fermi 
gases has attracted a great deal of interest recently [36l - l38| . This paper considers the three-body analog within a 
microscopic quantum mechanical framework. We note that our framework readily generalizes to bosonic three-body 
systems. The study of the dimensional crossover of bosonic systems is interesting as it allows one to study how, under 
experimentally realizable conditions, Efimov trimers (soj that are known to exist in three-dimensional space disappear 
as the confinement geometry is tuned to an effectively low-dimensional geometry [ioj . 

This paper generalizes the methods developed in Refs. 8, 32] for three equal-mass fermions in two different pseu- 
dospin states under spherically symmetric harmonic confinement to anisotropic harmonic confinement. We develop 
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an efficient and liiglily accurate algoritfim to calculate the eigenenergies and eigenstates of the system up to relatively 
high energies as functions of the interaction strength and aspect ratio of the trap. Several applications are considered: 
(i) The BCS-BEC crossover curve is analyzed throughout the dimensional crossover, (ii) For large and small aspect 
ratios, the energy spectra are analyzed in terms of strictly one-dimensional and strictly two-dimensional effective 
three-body Hamiltonian. (iii) The second- and third-order virial coefficients are analyzed as functions of the temper- 
ature, aspect ratio and scattering length. In particular, we show that the high-temperature limit of the third-order 
virial coefficient 63 at unitarity is independent of the shape of the trap in agreement with expectations derived through 
use of the local density approximation. For finite scattering lengths, 62 and 63 for anisotropic harmonic confinement 
are well approximated by those for isotropic harmonic confinement. 

The remainder of this paper is organized as follows. Section [II] presents a formal solution to the problem of three 
s-wave interacting fermions confined in an axially symmetric harmonic trap. We also consider the extreme cases of 
strictly one-dimensional and strictly two-dimensional confinement. Sections IIIII and IIVI applv the formal solution to 
cigar-shaped and pancake-shaped traps, respectively. We determine a large portion of the eigenspectrum as a function 
of the scattering length and discuss the transition to strictly one-dimensional and strictly two-dimensional geometries. 
Section |V] uses the two- and three-body eigenspectra to calculate the second- and third-order virial coefficients as a 
function of the temperature and the geometry of the confinement. Finally, Sec. IVII concludes. 



II. FORMAL SOLUTION 



We consider a two-component Fermi gas consisting of two spin-up atoms and one spin-down atom with interspecies 
s-wave interactions under anisotropic harmonic confinement. We refer to the two spin- up atoms as particles 1 and 2, 
and to the spin-down atom as particle 3. We introduce the single-particle Hamiltonian HQ{rj,A4) for the j*'' particle 
with mass A4 under harmonic confinement, 

Hoir,,M) = ^Vl^ + \m{ujIz^^+loIp]). (1) 

Here, is measured with respect to the trap center, and in cylindrical coordinates we have = (zj, pj, 0^). In 
Eq. ([T]), LOz and ujp are the angular trapping frequencies in the z- and p-directions, respectively. The aspect ratio 
rj of the trap is defined through rj = ojp/ujz- In this paper, we consider cigar-shaped traps with > 1 as well as 
pancake-shaped traps with 77 < 1. Our three-particle Hamiltonian H then reads 

3 

H = Y,Ho{vj,M) + Vint, (2) 

where Vint accounts for the interspecies s-wave two-body interactions, 

M„t = V^pf (r3i)-f l^pf (r32). (3) 
The regularized pseudopotential V^^ is characterized by the three-dimensional s-wave scattering length [23 - [26j , 



where Vjk = - ffc and r^k = |rjfe|- 

Since the trapping potential is quadratic, the relative and center of mass degrees of freedom separate and we rewrite 
the Hamiltonian H in terms of the relative Hamiltonian H^^i and the center of mass Hamiltonian Hem, H — Hici+ Hem- 
In the following, we obtain solutions to the relative three-body Schrodinger equation i/rci^' — -£'3b^, where 

Hrcl = Hjcifi + Vint (5) 

with 

H,cio = Ho{r, + HoCR, /x). (6) 

In Eq. ([5]), ^ is the two-body reduced mass, ^ — Ai/2, and the relative Jacobi coordinates r and R are defined 
through r = r^i and R — ^ C"^^*"^ — r2). Depending on the context, we use either r and R or r3i and r32 to describe 
the relative degrees of freedom of the three-body system. 
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To determine the relative three-body wave function ^'(r,R), we take advantage of the fact that the solutions to 
the "unperturbed" relative Hamiltonian i?ici o ^re known and consider the Lippmann-Schwinger equation (see, e.g., 
Ref. 0) 

*(r, R) - - / G(S3b; r, R; r', R')Mnt(r', R')*(r', R') dr'dTL'. (7) 

The Green's function G for the two "pseudoparticles" of mass fi associated with the Jacobi vectors r and R is defined 
in terms of the eigenstates $Ai (r)<I'A2 (R) and the eigenenergies Ex^ + of i?rcLo, 

$* (r')$t (R')$;^,(r)$A,(R) 
G(^3b;r,R;r^R0 = y ■ (8) 

Here, A collectively denotes the quantum numbers needed to label the single-particle harmonic osillator states. In 
cylindrical coordinates, we have A — (ri^, rip, m) with = 0, 1, 2, • • • , np = 0, 1, 2, • • • , and to = 0, ±1, ±2, • • • . The 
single-particle harmonic oscillator eigenenergies and eigenstates read 

Ex = (n^ + l]haj^ + {2np + \m\ + 1) tjHuj^ (9) 



and 



where 



and 
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$A(r)-^„,(z)i?„^,™(p)-=, (10) 



= \/a../^2- n -^"P 1-2^ (11) 



- J^IT^^ -p f . (12) 



alin, + \m\)l ^ V 

In the last two equations, Hn^{z/az) and Ln™'^ (rjp^ ja?^ denote Hermite and associated Laguerre polynomials, re- 
spectively. Throughout most of Sees. UlllIVl we use the oscillator energy Ez and oscillator length az [Ez = fiojz and 
Uz = \/hl (/iWz)] as our energy and length units. 

In Eqs. ([71)- (|12l) . we employ cylindrical coordinates since this choice allows us to write the Green's function G 
compactly. However, the two-body s-wave interaction potential is most conveniently expressed in spherical coordinates 
[see Eq. Since the pseudopotential V^^g°(r) acts only at a single point, namely at r = 0, it imposes a boundary 
condition on the relative three-body wave function \I'(r,R) (see, e.g., Ref. [3J]), 



*(r,R)Uo-^(^-4^)- (13) 



The unknown function /(R) can be interpreted as the relative wave function of the center of mass of the interacting 
pair and the third particle. Similarly, the pseudopotential V^^{v2,-i) imposes a boundary condition on the wave 
function \l/(r,R) when 0. Since the wave function \['(r,R) must be anti-symmetric under the exchange of 

the two identical fermions, i.e., Pi2^'(r, R) = — ^(r, R), where P\2 exchanges particles 1 and 2, the properly anti- 
symmetrized boundary condition corresponding to Vp^s"'^(r32) reads 



V&(r32R32)l ^_ /(R32) (^_^ 



(14) 



Here, we defined R32 = ^ (^^^ - ri) . 

To simplify the right hand side of Eq. (O, we impose the limiting behaviors of ^(r',R') for r'^-^ — ^ and — > 0, 

and expand /(R') in terms of the non-interacting harmonic oscillator functions, /(R') — X^a' /^''J'a' (R')- Using 
Eq. ([5]) for G and orthonormality of the single-particle harmonic oscillator functions, we find 



3/2 

EmJ 



X 

X 



Q'"" (i?3b - i?A; r; 0) $a(R) - g^"" i i?3b - Ey, '-^^M $a ( ^^^^ 



(15) 
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Here, we used that can be written as (r + -\/3R-)/2 and introduced the one-body Green's function (i?; r; r') 
for the pseudoparticle of mass /i that is associated with the relative distance vector r' , 

The one-body Green's function Q^^ (£'2b;r;r') with r' = coincides with the solution to the relative Schrodinger 
equation for two particles under harmonic confinement interacting through the zero-range pseudopotential Vp^{r) 
with s-wave scattering length a^^ and relative two-body energy E2h- {E; r; 0) is known for all aspect ratios ?? 
(see also Sees. Illll and IIV|) . 

To determine the expansion coefficients fx, we apply the operation ^(f-)|r-j.o ^^^^ hand side and the right 

hand side of Eq. (|T5|) . i.e., we multiply both sides of Eq. p5|) by r, then apply the derivative operator and lastly take 
the limit r — !■ 0. Defining 



J-3D(eA, ?7) = 2Tr E,al ^ {r^^^ ([ex + V+ 1/2]^^.; r; 0)} 



(17) 



with {ex + ri + l/2)E^_ = E31, - Ex, we find 

-^E/v*A'(R) = 



27ra3D 

A' 



E fy I ^•^'°(^V, ^)$A' (R) - E.alG'^ ^[ey + ri + 1/2]E.; ^R; 0^ $v 



(18) 



The quantity e^' can be interpreted as a non-integer quantum number associated with the interacting pair. If we 
multiply Eq. (I18p by $^(R) and integrate over R, we find an implicit eigenequation for the relative three-body energy 
i?3b or equivalently, the non-integer quantum number ex, 



n 

a' 



where 



iT^x' (^A') = 2nE^al J G^^ \^[ey + + 1/2]E,; ^R; j $x' $1(R) dK (20) 

and Sxx' is the Kronecker delta symbol. The determination of ^x^'(eA') ^-^^d J^^^i^x' tV) for 77 > 1 and 77 < 1 is 
discussed in Sees. IIIII and IIVI respectively. 

Equation (|19p can be interpreted as a matrix equation with eigenvalues a^/a'^^ and eigenvectors fx [1, [l3]- In 
practice, we first calculate the matrix elements ^x^'(ex') ^l- (HH) fo^' ^ given three-body energy i?3b and obtain the 
corresponding scattering lengths for this energy by diagonalizing the matrix with elements ^^'^'(ex')~-^^'^(^A, v)^x,x' ■ 
This step is repeated for several three-body energies. Lastly, we invert a'^°(£'3b) to get i?3b(a'^°), i.e., to get the three- 
body energies as a function of the s-wave scattering length. 

Equation (|19p has a simple physical interpretation. If the interaction between particles 2 and 3 is turned off, the 
matrix -^x'a' (^-'^' ) vanishes and the solution reduces to that of an interacting pair (particles 1 and 3) and a non- 
interacting spectator particle (particle 2). The relative energy (ex' + r] + l/2)Ez of the pair is determined by solving 
the relative two-body eigenequation J^^^{ex' ,11) — —Uz/a^^. The matrix -^x^/(ex') th^s arises from the fact that 
particle 3 not only interacts with particle 1 but also with particle 2. Correspondingly, the terms in Eq. (I19p that 
contain ^x^'(^-'^') '^^^ interpreted as exchange terms that arise from exchanging particles 1 and 2 [l3 |. 

For ry = 1, the function J^^^{€x, v) is given in Table |T] and the evaluation of ^x^'(eA') been discussed in detail in 
Ref. [T5I. The r] ^ 1 cases are discussed in Sees. Illll and IIVI For a spherically symmetic system with rj = 1, the total 
relative angular momentum quantum number L, the corresponding projection quantum number M and the parity 
n are good quantum numbers, and the eigenvalue equation can be solved for each L and M combination separately 
using spherical coordinates [8|. For a fixed L and M, A = (n, I, m) and A' = (n', I' , m') in Eq. (IT^ are constrained by 
I ^ V ~ L and m = m! — M. The parity of the three-body system is given by 11 = (—1)^. 
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TABLE I: Two-body properties in three dimensions (s-wave channel), two dimensions (m = channel), and one dimension 
(even parity channel). ip{ci) denotes the relative two-body wave function and fi the two-body reduced mass, q stands for r, p, 
and 2 in three, two, and one dimensions, respectively. For 3D, 2D and ID, we have e — E^h/Ez — f) — 1/2, e = E2h/Ep — 1 and 
e = E2h/Ez — 1/2, respectively, where E2h denotes the relative two-body energy, i/jg denotes the digamma function and 7 the 
Euler constant, 7 ~ 0.577. 



3D 



2D 



ID 



?^°'5'^'(r)^r 



/"^'^^(P)|;P 



An — a 
n 



Inf 



. + 1]- 



V'(q)l 



oc [ln(a™) - ln(p)] 



oc (2 - 



Bethe-Peierls B.C. 



a(r-i/.) 
dr 



ln(p/a2D) 



J" 



J^^"(e) = TvEpC 



{t [rQ'"" {[^ + n+ \\E.; r; 0)] {[e + l]Ep; p; 0) + ln(p/a,)}^^^, {G^^ {[e + \]E.;z; O) }^ 



-F='°(.,l) = 



-2r(-6/2) 
r(~./2-i/2) 



.F^°(6) = -iv>,(-6/2)-7 



r(-<:/2) 
2r(-e/2 + l/2) 



two-body energy F {e,ri) = —Uz/a 



We emphasize that the outlined formalism makes no approximations, i.e., Eq. ()19p with i^x') given by Eq. ([20)1 
describes all eigenstates of ffici [see Eq. ([SJ] that are affected by the interactions. In particular, all "channel couplings" 
are accounted for. In practice, the construction of the matrix -^^^'(^a') requires one to choose a maximum for A and 
A', or alternatively, a cutoff for the single-particle energy Ex. As has been shown in Ref. [11], good convergence is 
achieved for a relatively small number of "basis functions" for 77 = 1. As we show below, good convergence is also 
obtained for anisotropic confinement geometries. 

The formalism outlined can also be applied to strictly one-dimensional and strictly two-dimensional systems. Ta- 
ble H] defines the one-dimensional and two-dimensional pseudopotentials as well as a number of key properties of the 
corresponding relative two-body system. Making the appropriate changes in the outlined derivation and using the 
properties listed in Table HI we find for strictly one-dimensional systems 



E 

n'=0 



=0 

where is defined in Table HI 

IZk M = ^.«./" (^[^K + imE.;^z;0j (-0 <(z)dz, (22) 

and i?3b — En' = (e„' -I- l/2)£'2. Here, i?„^ denotes the single-particle energy of the one-dimensional system, i?„^ — 
(n^ -h 1/2)^;^, and Q^^ {E; z; z') the one-dimensional even parity single-particle Green's function, 

(23) 

For z' = 0, the single-particle Green's function is given by 

(E: 0) = ,^»p r i-^J^^\ U (-^^. i i ) , (24, 



2^E^a^ ' \ 2alJ \ 2 J \ 2 ' 2' a2 

where T(x) is the Gamma function and C/(a, 6, z) the confluent hypergeometric function. The strictly one-dimensional 
relative three-body wave function ^E* is characterized by the parity Hz- For even parity states, i.e., for states with 
Hz = 1, Uz and n'z in Eq. (j2ip have to be even. For odd parity states, i.e., for states with Hz = —I, Uz and have 

to be odd. 

Similarly, for strictly two-dimensional systems, expressed in units of Ep and Op [Ep — fnup and Op = \/h/ (/iWp)], 
we find, in agreement with Ref. [16.] , 

CO 

\j-np,n'^,m{^n'^,rn) — ^ {^11^,111)^11^,11'^ /ri^ ,m = In ^ /rip ,m , (25) 

n' =0 
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where is defined in Table HI 

IZn'm{^n'n^)^{-lr^Ey^ I I [e,,, + l]E,; ^ p; \ R,.,^ [ ^ ] R,,^^„, (p) pdp, (26) 



J ^™ ( [en;,,m + l]Ep; —p; j i?„;^,„ (^-j Rn^,m (p) pdp, 



and -Bab — En'^^m = {^n'^.m + ^)Ep. Here, En^^m denotes the single-particle energy of the two-dimensional system, 
Enp,m — i'^'Up + |m| -I- l)Ep. The two-dimensional single-particle Green's function Q^^ {E; p; p') is defined analogously 
to the three- and one-dimensional counterparts [sec Eqs. ([T5| and (|23p ]. For p' = and states affected by the 
zero-range s-wave interactions 0] , one finds 

S-(.;«0) ^ ^o.p (-^) r r-S^k^) . (-^. u^) . (27, 



27rEpal \ 2alJ \ 2 J \ 2 ' ' a^y 

The strictly two-dimensional relative three-body wave function is characterized by the projection quantum number 
M and the parity Hp, Hp = (—1)^. For a fixed M, m in Eq. is constrained to the value m = M. The 
next two sections analyze, utilizing our results for strictly one- and two-dimensional systems, Eq. for cigar- and 
pancake-shaped traps. 

III. CIGAR-SHAPED TRAP 

To apply the formalism reviewed in Sec. |ll]to axially symmetric traps, we need the explicit forms of the functions 
G^^ ([ca + + 1/2] -Ez; r; 0) and J"'^^(eA, ??), that is, the relative solutions to the trapped two-body system. For cigar- 
shaped traps (77 > 1), it is convenient to write as H, 0| 

G'"" {[ex + V + mE.;r; 0) = ^ exp i^-^j {vpV<^l) ([^x - 2^ + l/2]i?.; z- 0) , (28) 

where G'^^ {E; z; 0) is defined in Eq. Using Eq. ^ in Eq. 1^, we obtain 

/3D,(6v) = v^(-l)'"<5™,™' X 

jmax 

lim V (eA',jVnp,n' (29) 

where 

Il,r.'My,j)^E,a, I G'"" I [ey-2r^j + l/2]E,;^z;0 ] (-1) ^„^(z) (30) 



j=0 



{ey,j)^E,a,J _2ry.? + l/2]i?z;^^;oj (-|) ^„,(z) 



and 



/'^ ' (7) 

rip ,m / 







02 / Rjfi I — P I i?n;,m j Rrip,m{p) pdp. (31) 



The evaluation of the integrals /^^^„/ (^a'' j) ^-^d /j^j^.n' discussed in Appendix lAl The superscript "c" indicates 

that the integrals apply to cigar-shaped systems; for pancake-shaped systems (see Sec. lIVp . we introduce the integrals 
^np.n;,m(ev,j) and /P^_„, (j) instead. 

Although it is possible to calculate J-^^ {e\,ri) numerically for any trap aspect ratio 77, we restrict ourselves to 
integer aspect ratios for simplicity. For traps with integer aspect ratio, an exact analytical expression for J-"^'^(eA, ^7) 
is known 1, 0] , 



E..rri,-f;i-^;expf^)), (32) 



k=l 
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FIG. 1: Relative three-body energies E^^^/Ez as a function of the inverse scattering length azjuC"^ for a cigar-shaped trap with 
aspect ratio rj — 2 and (a) AI — and fl^ = -f 1, and (b) AI — and Hz — —1. 

where 2-Fi(a, b; c; z) is the hypergeometric function [47*1. Knowing -^^^'(eA') ^^'^ -^^^i^x'j v)^ Eq. p9)) can be diagonal- 
ized separately for each {Hz, M,Ilp) combination. We recall from Sec. [TTlthat A — {nz,np,m) and A' = (n'^, nj,, m'). 
The m and m' values are constrained by m — m' = M . Moreover, for — +1 and XI^ = —1, we have = = even 
and riz — n'z — odd, respectively. 

Figure [1] shows the three-body relative energies E^i^/Ez for 77 = 2 for states with (a) M = and = 4-1 and (b) 
M — Q and 11^ = — 1 as a function of the inverse scattering length a^/a'^^. The non- interacting limit is approached 
when {a^^)^^ — > ±00, and the infinitely strongly-interacting regime for (a'^^)^^ = (center of the figure). For 
each fixed projection quantum number M, we include around 840 basis functions. This corresponds to a cutoff of 
around (82 -I- 2M)Ez for the single-particle energy Ex. We find that jmax ^ 30 yields converged values for ^^^'(eA')j 

Eq. (f29|) . For small la^^/ozl (a^° positive and negative), our eigenenergies agree with those obtained within first-order 
perturbation theory. Our analysis shows that the energy of the ground state at unitarity has a relative error of the 
order of 10~^. The accuracy decreases with increasing energy. For example, for energies around 20Ez, the relative 
accuracy at unitarity is of the order of 10~^. 

The eigenstates fall into one of two categories: atom-dimer states and atom-atom-atom states. The eigenenergies 
associated with the former are negative for large positive Qz/a^^ while those associated with the latter remain positive 
for large positive Uz/a^^. The energy spectra shown in Fig. [T] exhibit sequences of avoided crossings. To resolve these 
crossings, a fairly fine mesh in the three-body energy is needed. In the {a^^)^^ — ^ —00 limit, the lowest M = state 
has negative parity in z, i.e., 11^ — —1. This is a direct consequence of the fact that the two identical fermions cannot 
occupy the same single particle state. In the {a^^)~^ — ?► -l-oo limit, in contrast, the lowest M = state has positive 
parity in z, i.e., II^ = -1-1. This is a direct consequence of the fact that the system consists, effectively, of a dimer and 
an atom. 

The main part of Fig. [2] shows the relative energy of the energetically lowest-lying state, the so-called crossover 
curve, of the three-body system with M = for various aspect ratios of the trap (77 = 2, • • • , 10) as a function of 
the inverse scattering length az/a?^ . For comparative purposes, we subtract the ground state energy of 2r]Ez of the 
strictly two-dimensional non-interacting system, that is, the energy that the system would have in the p-direction if 
the dynamics in the tight confinement direction were frozen, from the full three-dimensional energy. In Fig.[51 asterisks 
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FIG. 2: (Color online) "Crossover curve" of the three-body system with M — 0, shifted by 2Ep — 2riEz, for cigar-shaped traps 
with rj = 2 (solid line), 4 (dotted line), 6 (dashed line), 8 (dash-dot-dotted line), and 10 (dash-dotted line) as a function of 
Oz/o^^. The scattering lengths at which the parity of the corresponding eigenstate changes from II^ = — 1 ("left side of the 
graph") to IIz = -1-1 ("right side of the graph") are marked by asterisks. At these points, the derivative of the crossover curve 
is discontinuous; the discontinuities are not visible on the scale shown. The inset shows the (unshifted) crossover curve as a 
function of ap/a^°. 



mark the scattering lengths at which the eigenstate associated with the crossover curve changes from Hz = — 1 to 
= -1-1. With increasing 77, the parity change occurs at larger a^/a^^ (that is, smaller a'^^/a^). The inset of Fig. [2] 
replots the crossover curves as a function of Up/a^^. 

We now discuss the large 77 limit in more detail. Using the limiting behavior of J^'^^(e, rj) for 77 3> 1 and rj ^ |e| [11,0, 



2,7J-^^(e) + V^C(l/2), 



the two-body eigenequation for the relative energy becomes 0, 01 
where the renormalized one-dimensional scattering length is given by [20l . [2]| 



"Ten 

a. 



,1D 



C(l/2) 



2a 



3D 



(33) 



(34) 



(35) 



Figure ^a.) shows the relative two-body energies for a system with 77 = 10, M = and Hz ~ +1 obtained by 
solving the eigenequation J^^°(e, 77 = 10) = — a^/a'^^ [see Eq. ([5^ for J^^^{e,ri)]. Figure (Sljb) compares the full 
three-dimensional energy (solid line) with the energy obtained by solving the strictly one-dimensional eigenequation, 
Eq. ([M)) . with renormalized one-dimensional scattering length aj^^ (dotted line). To facilitate the comparison, we add 
the energy of the tight confinement direction to the energy of the one-dimensional system. The agreement is quite 
good for all scattering lengths. The inset of Fig.[3l^b) shows the difference between the strictly one-dimensional energy 
and the full three-dimensional energy as a function of Uz/a^^. The maximum deviation occurs around unitarity and 
is of the order of 0.2%. 

Next, we discuss the behavior of the three-body system in the large 77 limit. If we use Eqs. (p3)) and (|35|) in Eq. (|19l) . 
we find 



E 

X' 



[2r, 



,1D 



' f 

-jx 



(36) 



A straightforward analysis shows that Eq. (I36|) reduces to its strictly one-dimensional analog if (i) the sum over A' 
is restricted to a sum over n'z [A' = (71^,0,0)]; (ii) the index A is restricted to A = (71^,0,0); (in) the energy E^y, 
is replaced by i?3b — 2i]Ez; and (iv) jmax in Eq. ([^ is set to zero. Under these assumptions, Eq. ([55)) reduces to 
Eq. (|2ip with a^^ replaced by a^^. We emphasize that the assumption 77 ^ \€x\ [see discussion around Eq. (p3)) ] is 



9 



W 15- 





' ... 












(a) . 

1 



-20 



12.5 



10.5 



a / a 

z 



3D 



20 



a / a 



3D 



20 











1 




.0.08 














1 








-0.04 
















A 






(b) - 


-loo 





loq 



20 



FIG. 3: (Color online) (a) Relative two-body energies E2h/Ez as a function of the inverse scattering length a^/a^^ for a cigar- 
shaped trap with aspect ratio i] = 10, M — and — +1. (b) The solid curve from panel (a) is replotted and compared with 
the energy obtained by solving the strictly one-dimensional eigenequation with renormalized one-dimensional scattering length 
OriJi (dotted line; the energy Ep = rjE^ has been added to allow for a comparison with the full three-dimensional energy). The 
inset shows the difference between the dotted and solid lines as a function of the inverse scattering length a^/a^^ . The scale 
of the y-axis is identical to that of the inset of Fig. \M)^)- 



not valid when two atoms form a tight molecule. In this limit, the three-dimensional s-wave scattering length, or the 
size of the dimer, is smaller than the harmonic oscillator length in the transverse direction, which implies that the 
strictly one-dimensional description is not valid. 

Figure lll^a) shows the relative three-body energies for states with M = and II^ = -t-1 as a function of the inverse 
scattering length a^/a^^ for a cigar-shaped trap with rj = 10. Figure |3lb) compares the energy of the energetically 
lowest-lying three-atom state with 11^ = -1-1 (solid line) [see thick solid line in Fig. |H[a)] with the corresponding 
state obtained by solving the strictly one-dimensional equation with renormalized one-dimensional scattering length 
aj^j. We also include the energy of one of the eigenstates with M = and 11^ = —1 (dashed hne). The inset 
shows the difference between the energy obtained within the strictly one-dimensional and the full three-dimensional 
frameworks. The maximum of the deviation occurs near unitarity. The agreement between the full three-dimensional 
and the strictly one-dimensional descriptions is good. Importantly, the deviations for the three-body system with 
M = and Hz — +1 [solid line in the inset of Fig. lUJb)] are only slightly larger than those for the two-body system 
[inset of Fig. [3l^b)], suggesting that the presence of the third atom does not, in a significant manner, reduce the 
applicability of the strictly one-dimensional framework — at least for states in the low-energy regime characterized as 
gas-like three-atom states. 
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FIG. 4: (Color online) (a) Relative three-body energies E^h/Ez as a function of the inverse scattering length a^/a'^^ for a 
cigar-shaped trap with aspect ratio = 10, M = and Hz — +1. (b) The solid curve from panel (a) is replotted and compared 
with the energy obtained by solving the strictly one-dimensional eigenequation with renormalized one-dimensional scattering 
length (the energy of 2Ep = 2t]Ez has been added to allow for a comparison with the full three-dimensional energy). For 
comparison, the dashed line shows one of the three-dimensional energy curves for M = and Hz = — 1 [not shown in panel (a)]; 
the corresponding strictly one-dimensional energy is shown by a dotted line. The difference between the full three-dimensional 
and strictly one-dimensional descriptions is hardly visible on the scale shown. Solid and dashed lines in the inset show the 
differences between the strictly one-dimensional energies [dotted lines in panel (b)] and the full three-dimensional energies [solid 
and dashed lines in panel (b)] as a function of the inverse scattering length a^^joF^ for Hz = -fl and Hz — —1, respectively. 



IV. PANCAKE-SHAPED TRAP 



For pancake-shaped traps with 77 < 1, we use the following form of the Green's function [l, 0, 



^-([e. + . + l/2]i..; r; 0) = ^ exp (-^) g (./a.)^- 



(37) 



This expression is equivalent to Eq. psp but converges faster for pancake-shaped traps than Eq. psp . Using Eq. ([37l 
in Eq. (1^ . we obtain 



/3D,(ev) = 2V^(-l)"<5„ 



j.^ tW^MIljP (,VP (e, ,) 



where 



1/2 



/V3 \ 



(38) 



(39) 
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FIG. 5: Relative three-body energies Eyo/ Ep as a function of the inverse scattering length ap/a?^ for a pancake-shaped trap 
with aspect ratio 77 = 1/2 and (a) M — and Hz = +1, and (b) AI — ±1 and II^ = +1. 



and 



00 

2 / ^2D 



Ezai / g 



2j 



V3 



(40) 



Details regarding the evaluation of the integrals are explained in Appendix [Xj In the following, we limit ourselves to 
cases where the reciprocal of the aspect ratio is an integer. In this case, we have [1, 0] 



E 



ex + l 



/cry) 



(41) 



Figure [5] shows the relative three-body energies E^\y/Ep as a function of the inverse scattering length for 77 — 1/2, 
+1, and (a) M = and (b) M = ±1. In the {a?^)~^ -00 hmit, the ground state has M = ±1 and 11^ = +1 
symmetry. In the {a^^)~^ — ^ -l-oo limit, in contrast, the ground state has M = and IIj = -1-1 symmetry. 

Figure [6] shows the relative energy of the energetically lowest-lying state, the so-called crossover curve, of the three- 
body system with II^ = +1 for various aspect ratios of the trap (77 — 1/2, • • • , 1/10) as a function of the inverse 
scattering length ap/a?^ . For comparative purposes, we subtract the relative ground state energy of Ez of the non- 
interacting one-dimensional system, that is, the energy that the system would have in the z-direction if the dynamics 
in the tight confinement direction were frozen, from the full three-dimensional energy. The scattering lengths at 
which the symmetry of the corresponding eigenstate changes from M — ±1 to M = are marked by asterisks. The 
symmetry change occurs around az/a?'^ ~ 1 (see inset). 

It is instructive to compare Fig. [5] (pancake-shaped trap) and Fig. [5] (cigar-shaped trap). For both geometries, the 
crossover curve changes symmetry. The change of the symmetry is associated with the low-energy coordinate (the 
p-coordinate for pancake-shaped systems and the z-coordinate for cigar-shaped systems). For both geometries, the 
symmetry change occurs, for the aspect ratios considered, when a^^ is of the order of the oscillator length in the tight 
confinement direction. 
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FIG. 6: (Color online) "Crossover curve" of the three-body system with XI^ = +1, shifted by — Ep/rj, for various aspect 
ratios of the trap, 77 = 1/2 (solid line), 1/4 (dotted line), 1/6 (dashed line), 1/8 (dash-dot-dotted line), and 1/10 (dash-dotted 
line) as a function of Up/a^^ . The scattering lengths at which the M quantum number of the corresponding eigenstate changes 
from M = ±1 ("left side of the graph") to M = ("right side of the graph") are marked by asterisks. The inset shows the 
(unshifted) crossover curve as a function of a^/a'^^ . 

Next, we consider the small ry limit in more detail. For rj <^ 1 and |e| <^ 1, we have [sl, 01 

« [2-^'°(^) - 21n(C) - ln(^)] , (42) 
and the two-body eigenequation for the relative energy becomes 

^'^ie)^Ha^Z/ap). (43) 
where the renormalized two-dimensional scattering length is given by [2^ 

^^y^Cexpf-^) (44) 



with C « 1.479. Figure [TJa) shows the relative two-body energies for a system with rj = 1/10, M — and = 
+1 obtained by solving the eigenequation J^'^^(e,7] ~ 1/10) = — a^/a'^^ [see Eq. (PT|) for J^^^ [e,!])]. Figure [71[b) 
compares the full three-dimensional energy (solid line) with the energy obtained by solving the strictly two-dimensional 
eigenequation, Eq. (j43l) . with renormalized two-dimensional scattering length (dotted line). For comparative 
purposes, we add the energy of the tight confinement direction to the energy of the strictly two-dimensional system. 
The inset of Fig. shows the difference between the strictly two-dimensional energy and the full three-dimensional 
energy as a function of ap/a^^ . The maximum deviation occurs around unitarity and is of the order of 0.4%. 
To treat the three-body system in the small rj limit, we insert Eqs. and (|44l) into Eq. ([T^ . This yields 



E 

A' 



^rv(^A')--^^"(ev)<5A.A' 



a 



h' - In ( Ti^ ) fx. (45) 

rcn 



For fixed M, Eq. P5|) reduces to the strictly two-dimensional eigenequation, Eq. (^5]) . if (i) the sum over A' is restricted 
to a sum over n'^ [i.e., if A' = {0,n'p,m' = M)]; (ii) the index A is restricted to A = (0, rip, to = M); (in) the energy 
i?3b is replaced by i?3b — E^; and (iv) jmax in Eq. ([551) is set to zero. Under these assumptions, Eq. (^51) reduces to 
Eq. ([25]) with 0^° replaced by 



2D 



Figure mja) shows the relative three-body energies for states with M = and ~ +1 as a function of the inverse 
scattering length Op/ a^^ for a pancake-shaped trap with r; = 1/10. Figure[8Ub) compares the energy of the energetically 
lowest-lying three-atom state with M = (solid line) [see thick solid line in Fig. HJa)] with the corresponding state 
obtained by solving the strictly two-dimensional equation with renormalized two-dimensional scattering length 
(dotted line). The inset shows the difference between the energies obtained within the strictly two-dimensional and 
the full three-dimensional frameworks. Similar to the one-dimensional case, the maximum of the deviation occurs 
near unitarity. Comparison of the insets of Figs.jTJb) and[5l^b) suggests that, at least in this low-energy example, the 
presence of the third atom does not, in a significant manner, reduce the applicability of the strictly two-dimensional 
framework. For the same aspect ratio, the deviations are expected to increase with increasing energy. 



13 



15 




...! ' 


Q. 






PJ 






IN 












5 






. w- 



-50 50 




FIG. 7: (Color online) (a) Relative two-body energies E2h/Ep as a function of the inverse scattering length Up/a^^ for a 
pancake-shaped trap with aspect ratio 77 = 1/10, M = and Hz = +1- (b) The solid curve from panel (a) is replotted and 
compared with the energy obtained by solving the strictly two-dimensional eigenequation with renormalized two-dimensional 
scattering length a^^^ (dotted line; the energy -Ez/2 has been added to allow for a comparison with the full three-dimensional 
energy). The inset shows the difference between the dotted and solid lines as a function of the inverse scattering length Up/a^^ . 
The scale of the j/-axis is identical to that of the inset of Fig. [Sjb). 



V. SECOND- AND THIRD-ORDER VIRIAL COEFFICIENTS 

This section utilizes the two- and three-body energy spectra to determine the second- and third-order virial co- 
efficients as functions of the s-wave scattering length a?^ , aspect ratio 77 and temperature T. The n*''-order virial 
coefficient 6„ enters into the high-temperature expansion of the grand-canonical thermodynamic potential VL of the 
equal-mass two-component Fermi gas with interspecies s-wave interactions where 
ri^i) and ri^^) denote the grand-canonical thermodynamic potential of the spin-up component and the spin-down 
component, respectively, and fi^i^-* accounts for the interspecies interactions, 

00 

n=2 

Here, z is the fugacity, z = exp[/i/(fcBT')]. In the high-temperature limit, z is a small parameter and the expansion 
given in Eq. (j46l) is expected to provide a good description if the sum is terminated at quadratic or cubic order. 
The coefficient 62 of the term is determined by one- and two-body physics and the coefficient 63 of the z"^ term 
is determined by one-, two- and three-body physics. As the temperature approaches the transition temperature 
from above, the de Broglie wave length increases and, correspondingly, the fugacity z increases. It follows that the 
number of &„ coefficients needed to accurately describe the thermodynamics increases with decreasing temperature. 
Comparison with experimental data has shown that the inclusion of &2 and 63 yields quite accurate descri ptio ns of 
the high-temperature thermodynamics of s-wave interacting two-component Fermi gases (see, e.g., Refs. Bllilll)- 
In Eq. (j46l) . Qifl denotes the canonical partition function of a single spin-up particle. We define the canonical 



14 




FIG. 8: (Color online) (a) Relative three-body energies Esh/Ep as a function of the inverse scattering length ap/a?^ for a 
pancake-shaped trap with aspect ratio r] = 1/10, M = and Hz = -1-1. (b) The solid curve from panel (a) is replotted and 
compared with the energy obtained by solving the strictly two-dimensional eigenequation with renormalized two-dimensional 
scattering length (dotted line; the energy of has been added to allow for a comparison with the full three-dimensional 
energy). The inset shows the difference between the dotted and solid lines as a function of the inverse scattering length ap/a?^ . 



partition function Qm,n2 of the system consisting of ni spin- up particles and 712 spin-down particles through 

Q„„„, = ^exp ' , (47) 

where (711,7x2) denotes the total energy of the system (including the center-of-mass energy) and the summation 
over j includes all quantum numbers allowed by the symmetry of the system. For equal-mass fermions, as considered 
throughout this paper, we have (5i,o = Qo,i = Qi- The virial coefficients 62 and 63 can be expressed as [13, El 

= 9l±_^ (48) 

and 

, r, Q2,i — Q2,oQi — b2Qi 



Qi 



(49) 



The virial coefficients 62 and 63 depend on the interspecies scattering length a^°, aspect ratio 77 and temperature 
T. Once the thermodynamic potential is known, physical observables such as the pressure and the entropy can be 
calculated. 

We now discuss the determination of 62 and 63 for equal-mass two-component Fermi gases under anisotropic 
harmonic confinement. In this case, the single-particle canonical partition function Qi can be determined analytically. 



exp([l/2-K77] Cj^) 
[exp (0)2) - 1] [exp (77^)2) - 1] 



^ . ^^^UV^-r-/J^7 .2 ^ (50) 
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where Cj^ denotes the "inverse temperature" in units of Ez^ Cjz = Ez/{k'QT). Alternatively (see below), we express 
the inverse temperature in units of Ep or Save, — Ep/{k-QT) and Wave(''7) = £'ave/(^BT')- The average energy i^ave 
is defined in terms of the root-mean-square or, in short, average angular frequency WaveC??), -^avc — ^avc(?7), where 



Wavo(?7) = Y ■ (51) 

We note that the average angular frequency coincides with the angular trapping frequency for rj = 1 but not for t] ^ 1. 
Below, we frequently suppress the explicit dependence of Wave on -q. The partition functions Qi^i and Q2.1 can be 
determined from the two- and three-body energy spectra (see Sees. IIIII and lIV[) for each s-wave scattering length a'^^, 
aspect ratio t] and temperature T. 

At unitarity, the high-temperature expansion of 6„ reads (for rj — 1, see Ref. [l3 |) 

K « bf^ + bi? + 6(f) {7^)C:t + ■■■. (52) 

The coefficients b''n \ that is, the high-temperature limits of the trapped virial coefficients bn, are independent of the 
aspect ratio rj. This has previously been shown to be the case for n = 2 Here, we extend the argument to all 
n. Application of the local density approximation [l^. [l9l. 143^ to axially symmetric confinement potentials shows that 
the virial coefficients 6j^°™ of the homogeneous system, which have been shown to be temperature-independent at 
unitarity plI - fTsj . are related to bn'^ through 



L horn 



i^/^foW. (53) 



Since Eq. ([53)1 holds for all 77, bn^ has to be independent of 77 for all n. 

The expansion coefficients b'^rf', fc = 2,4, • • • , parametrize "trap corrections", that is, corrections that arise due to 
the fact that the harmonic confinement defines a meaningful (finite) length scale. In fact, for 77 7^ 1, the confinement 

defines two length scales, suggesting that the bn \ fc = 2, 4, • • • , depend on 77. Equation ([5^ expresses the temperature 
dependence of 6„ in terms of the inverse temperature associated with the z-direction, regardless of whether 77 is greater 
or smaller than 1. Interestingly, it was shown in Ref. 42| that the dependence of b^\ fc = 2,4, • • • , on the aspect 



ratio can be parametrized, to a good approximation, in terms of the average inverse temperature a;avo(77), 

bf\i^)Co1^b^^\l)[u,Ui^)]K (54) 

Equation (j54p implies that the trap corrections for two-body systems with rj ^ 1 can be parametrized in terms of the 
trap corrections for the spherically symmetric system if the inverse temperature is expressed in terms of the average 
trapping frequency that characterizes the anisotropic system. 

We now illustrate that Eq. ([51)1 applies not only to 62 but also to 63. Figures [5K a) andlSJb) show the third-order 
virial coefficient at unitarity for systems with > 1 and < 1, respectively. The virial coefficients are plotted as 
a function of the inverse temperature expressed in units of the weak confinement direction, i.e., in terms of ui^ for 
77 > 1 and in terms of ujp for 77 < 1. In the high-temperature limit, 63 approaches a constant, confirming that 63" •* 
is independent of 77. The dotted lines show the third-order virial coefficient for 77 = 1, calculated using the average 
trapping frequency characteristic for the respective anisotropic system. Figure [9] illustrates that the third-order virial 
coefficient for anisotropic traps is approximated well by that for 77 = 1 with appropriately scaled angular frequency. 
The insets of Fig. [5] show that the third-order virial coefficients of the anisotropic system collape, to a very good 
approximation, to a universal curve over a surprisingly large temperature regime, i.e., down to temperatures around 
k^T « i?ave/2. We conjecture that bn \'>])'^z can be approximated quite well by bit\^)[u}a.ve{'r])]'' for n = 4, 5, • • • as 
well, as long as fc is not too large, i.e, as long as the temperature is not too low. 

Next, we discuss the behavior of 62 for finite s-wave scattering lengths. Figure [TO] shows a surface plot of 62 for 
77 = 1 as a function of the inverse scattering length Cave/a^^ (a^^ < 0) and the inverse temperature Wavc- Here, a^vc 
denotes the oscillator length associated with the average trapping frequency, Oavo = \/^/ il^'-^mc)- At unitarity, 62 
is only weakly- dependent on the temperature and approximately equal to 1/2 (see discussion above). The smallest 
inverse temperature Wavo considered in Fig. [TO] is 0.0003. For this inverse temperature, &2 is fairly close to 1/2 for all 
dave/o?^ shown. Thus, Fig. [TO] shows that the high-temperature limit of 62 is nearly independent of the scattering 
length. This behavior can, as we now show, be understood from the two-body energy spectrum. 

Figure [TlT a) shows selected relative two-body energies as a function of Oavc/a'^^ for the trapped system with 77 = 1. 
In the low-energy regime (solid line), the two-body energy changes by nearly 2£'avc for the scattering length range 
shown. As the energy increases (dashed and dotted lines show energies around 200i?avo and 20000£'avo, respectively). 
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FIG. 9: (Color online) (a) Third-order virial coefficient 63 at unitarity as a function of the inverse temperature iDz for 77 = 1 
(solid line), 2 (dashed line), and 3 (dash-dotted line), (b) Third-order virial coefficient 63 at unitarity as a function of the 
inverse temperature cDp for rj = 1 (solid line), 1/2 (dashed line), and 1/3 (dash-dotted line). For 7^ 1, 63 terminates at the 
inverse temperature of about 0.25 since our calculations include a finite number of three-body energies; obtaining the behavior 
of 63 in the high-temperature limit requires the inclusion of infinitely many three-body energies. For 77 7^ 1, dotted lines show 
63 for ri — 1, calculated using the average frequency of the respective anisotropic system. This approximate description is quite 
good. The insets of panels (a) and (b) show the same data as the main figure, but now as a function of a)avc as opposed to uiz 
and UJp. The insets show that the third-order virial coefficients for different 77 collapse to a universal curve for all rj (deviations 
arise in the low-temperature regime, i.e., for ojavc ^ 1). 




FIG. 10: Second-order virial coefficient &2 for the two-body system under isotropic confinement as a function of the inverse 
scattering length aavo/a'"^ (a^'^ negative) and the inverse temperature tiavc. The smallest uj^ve considered is 0.0003. In the 
liavc — >■ limit, 62 approaches 1/2 for all aavc/a'"^ (a''° < 0; see text for further discussion). 
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FIG. 11: (Color online) (a) Relative two-body energies E2h/Eip,a, shifted by 2(n — 1), for isotropic confinement {jj = 1 and 
s-wave channel) as a function of the inverse scattering length a^-^c/a?^ . (b) Relative two-body energies E2h/Ep, shifted by 
2(n — 1), for pancake-shaped confinement (77 = 1/10, A/ = 0, and li^ = -1-1) as a function of the inverse scattering length 
ap/a^^ . Solid, dotted and dashed lines show the energies for n — 1, 100 and 10000, respectively. 



the two-body energy undergoes less of a change and eventually becomes nearly flat over the scattering length region 
shown. This implies that the high-energy portion of the two-body spectrum looks like that of the unitary gas over 
an increasingly large region around unitarity. The behavior of the energy spectrum can be understood by expanding 
the transcendental two-body eigenequation J-^^{e,l) = —aavc/ci^^ around unitarity. Assuming |aavc/a^^| 1, we 
find 4^ iij 

^2.+ n«.-^>±l/^^. (55) 



£^avc V 2/ nr{n + 1) 

Equation (1551) provides a good description of the energies as long as the right hand side is small. Since the right 
hand side of Eq. ([55]) scales for large n as {aave/cL^^)/\/n, Eq. ([55)1 provides, as n increases, a good description for 
an increasingly large region around unitarity. That is, the energies vary approximately linearly with Qavc/o-^^ , with a 
slope that approaches zero, as n — 00. This analysis rationalizes why 62 approaches 1/2 as T 00, regardless of the 
value of the scattering length (a"^^ < and finite). 

Figure [TU] has been obtained for a spherically symmetric system, that is, for rj — 1. We now demonstrate that 
Fig. [To] applies, for experimentally relevant temperatures, to all aspect ratios and not just to 77 = 1. Figure [T2l shows 
62 as a function of Wavc for three different aspect ratios, 77 = 1 (solid lines), ry = 1/5 (dashed lines) and r/ = 1/100 
(dotted lines). Three different scattering lengths are considered: a^^ = — Oave, 00 and aavc. It can be seen that 62 
is, to a very good approximation, independent of rj in the high-temperature (small Wavc) hmit. We stress that the 
independence of 62 of rj requires that the inverse temperature and scattering length are expressed in terms of the 
average oscillator units i?avo and a^vc, respectively. 

To understand the universality implied by Fig. [T^l that is, the fact that Fig. [TU] applies to all aspect ratios and not 
just to ?7 = 1, we analyze the behavior of the high- lying part of the relative two-body spectra for ?7 7^ 1. Figure [TTT b) 
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FIG. 12: (Color online) Second-order virial coefficient 62 of the trapped two-body system as a function of the inverse temperature 
liavc for three different values of the inverse scattering length (aavc/a"^° = — 1, and +1; see labels) for rj — I (solid lines), 
77 = 1/5 (dashed lines) and rj = 1/100 (dotted lines). 



exemplarily shows the relative energies for a pancake-shaped system with r; = 1/10. Comparison with Fig. Illf a) 
shows that the quahtative behavior of the high-energy part of the spectrum is independent of 77. This is confirmed 
by a more quantitative analysis that Taylor expands the implicit eigenequation for the anisotropic two-body system 
around unitarity. We conclude that two two-body systems with different aspect ratios but identical Oave/a^'^ and 
'^avc are characterized by approximately the same &2 value. Our analysis of the three-body energies for anisotropic 
confinement suggests that analogous conclusions hold for 63. We speculate that the conclusions hold also for the virial 
coefficients with n > 3. 

To estimate the extent of the universal behavior, it is instructive to reexpress uJavc in terms of the Fermi tempera- 
ture. For a spin-balanced system of N fermions under spherically symmetric confinement, we use the semi-classical 
expression fcBTp = {3Ny/^hw^^^, yielding Wavo = (37V)-i/3(t/Tf)-^ For N = 10^, 10^ and 10^, T/Tp = 1 cor- 
responds to (Z^avo ~ 0.149, 0.032 and 0.007, respectively. For these temperatures, the thermodynamic behavior is, 
according to our discussion above, expected to be essentially universal over a fairly wide range of scattering lengths. 
For T/Tp — 1, we estimate that the deviation of 62 from the value 1/2 approaches 5% for aave/a"^^ ~ —0.16, —0.38 
and —0.78 for N = 10^, 10"* and 10^, respectively. This implies that uncertainties of the scattering length dependence 
on the magnetic field in recent experiments at unitarity [13, [13 should have a negligibly small effect on the equation 
of state at unitarity. For all three cases considered above, the corresponding (/cpa'^^)"^ values is approximately —0.03. 

The determination of the high-temperature behavior of the third-order virial coefficient for different scattering 
lengths and aspect ratios is much more demanding than that of the second-order virial coefficient. Although our 
analysis of 63 is less exhaustive than that of 62, it suggests that the conclusions drawn above for the second-order 
virial coefficient carry over to the third-order virial coefficient. 



VI. CONCLUSION 



This paper developed a Lippmann-Schwinger equation based approach to determine the energy spectrum and 
corresponding eigenstates of three-body systems with zero-range s-wave interactions under harmonic confinement 
with different transverse and longitudinal angular trapping frequencies. The formalism was applied to the equal-mass 
system consisting of two identical fermions and a third distinguishable particle in a different spin-state. The energy 
spectra were determined as a function of the interspecies s-wave scattering length for various aspect ratios 77, ?/ > 1 
(cigar-shaped trap) and 77 < 1 (pancake-shaped trap). For ?/ 3> 1, we showed that the low-energy portion of the energy 
spectra are reproduced well by an effective one-dimensional Hamiltonian with renormalized one-dimensional coupling 
constant. Similarly, for 77 ^ 1, we showed that the low-energy portion of the energy spectra are reproduced well by an 
effective two-dimensional Hamiltonian with renormalized two-dimensional coupling constant. As the energy increases, 
the description based on these effective low-dimensional Hamiltonian deteriorates. 

The two- and three-body energy spectra were then used to determine the second- and third-order virial coefficients 
that determine the virial equation of state, applicable to two-component Fermi gases at temperatures above the Fermi 
temperature. Our key findings are: 
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(i) At unitarity, the second- and third-order virial coefficients 62 and 63 approach constants in the high-temperature 
regime. The constants (referred to as 62"'' and 63°^) are independent of ry. 

(a) For finite scattering length a^^ , we find that 62 and 63 cofiapse, to a very good approximation, to a single curve 
for all ry if the temperature and scattering length are expressed in terms of the average energy i^ave and the average 
oscillator length aavej respectively. Deviations from the universal curve arise in the low-temperature regime where the 
virial equation of state is not applicable. 

(in) The virial coefficient 62 is approximately equal to 1/2 over a fairly large temperature and scattering length 
regime around unitarity. 

The work presented in this paper is directly relevant to on- goin g cold atom experiments. The three-body spectra, 
e.g., can be measured experimentally through rf spectroscopy 2^ 4^|4^. Moreover, the formalism can be employed 



to characterize the molecular states in more detail, quantifying the "perturbation" of the dimer due to the third 
particle throughout the dimensional crossover. The determination of the virial coefficients is of immediate relevance 
to cold atom experiments that study the dynamics of large fermionic clouds under low-dimensional confinement. 
The formalism developed in Sees. HHITvI of this paper can be extended fairly straightforwardly to three-boson and 
unequal-mass systems. 
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Appendix A: Calculation of integrals involved 

This appendix presents the evaluation of the integrals defined in Sees. IIIII and HVl The integrals Inp,n' ,mU) ^^'^ 
n' U) energy-independent. They are tabulated once and then used for each three-body energy considered. To 
evaluate the integral /^^ „, „j(j), Eq. (pij) . we expand each of the three radial two-dimensional harmonic oscillator 

functions R [see Eq. ([T^ ] in terms of a finite sum, i.e., we use the series expansion of the associated Laguerre 
polynomials |47| . 

(n- j)!(fc + j)!j! 

7^ n' m(i) then becomes a finite triple sum that contains integrals of the form 

The finite sum is readily evaluated and /^^ „, is stored for each rip, n'p, m and j combination. To evaluate 

n' 0)' P^P - ■'^fi rewrite each of the three one-dimensional harmonic oscillator functions ip [see Eq. pTjl ] in 
terms of associated Laguerre polynomials instead of Hermite polynomials j48| , 

H2n{z/a,) = (-l)"22"n!4-i/2)(^2/^2) (A3) 

and 

i?2„+i(^K) = (-l)"22"+in!(zK)4i/2)(^2/^2)^ 

The evaluation of /^^ {j) then proceeds analogously to that of In^^n' .mU)- 

The integrals and mi^x'd)^ Eqs. (|30|) and (|40|) . are energy- dependent. To evaluate these inte- 

grals, we first "separate out" the energy dependence and then proceed along the lines discussed above. The energy 
dependence of 7,';^ „, (e^, j) enters through g^^{[ex' - 2rjj + l/2]E,; V3/2z; 0). Using the identity ^ 

00 ^(-1/2)/ N 

r(a)L/(a, 1/2, x) ^ ^ (A5) 

k=0 
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FIG. 13: Convergence study, (a) The dots show the coefficient C^^ , Eq. (IA7|l . for = 2 and n'^ = 80 as a function of k. 
(b) The circles and triangles show the sum JJ^^ Eq. (|A6|I . for — 2r]j)/2 — 25.745 and 40.234, respectively, as a 

function of the cutoff fcmax- As in panel (a), we used = 2 and n'^ = 80. 



^n- n' (^x'tJ) '^^n. be written as an infinite series, 

^ fcmax ^ 

where the coefficients „, are energy-independent, 

CLk = l\nA^)^Ki~z/2)exp (^-g^4-i/^) (^g^ d^. (A7) 

The evaluation of the C^^ „, now proceeds as above. 

The series introduced in Eq. (jA6P diverges if {e^' — ^vj)/'^ is a positive integer or zero, that is, in the non-interacting 
limit. In practice, this does not pose a constraint since the energy grid can be chosen such that it does not contain the 
non-interacting energies. To analyze the dependence of /^^ (e^'ii) tti^ upper summation index fcmax, Fig. fTSfa) 
shows the behavior of C^^ „, for fixed Uz and n'^ as a function of fc. In this example, C^^ „, vanishes nearly identically 
for k > 30. The coefficients C^^ are multiplied by the fc-dependent prefactor [fc — {ey — 2r]j)/2]~-^. This prefactor 
is maximal for fc « {ex' — 2'qj)/2. Figure \W(h) shows the quantity I^^^n' (^a'jJ) as a function of fcmax for the same 
Uz and n'z as in Fig. WSi&] but two different values of (e_^' — 2r]j)/2, that is, for {e^' ~ 2rij)/2 — 25.745 and 40.234. 
In the first case, the absolute value of the prefactor takes its maximum at fc = 26, where the coefficients C^^ „/ are 
non-zero. Correspondingly, /^^ shows a sharp peak near fcmax = 26 and then smoothly approaches its asymptotic 
value [see circles in Fig. [TBT b)]. In the second case, the absolute value of the prefactor takes its maximum at fc = 40, 
where the C^^ „, coefficients are very small. Correspondingly, the C^^ ^, coefficients suppress the maximum of the 
prefactor and the quantity /^^ approaches its asymptotic value for fcmax ~ 30 [see triangles in Fig. [TSl b)]. We find 
that the choice of fcmax ^ 2max(n2,n'2) yields well converged values for /^^ „, (ex', j) foi' all energies considered. 

In an alternative approach, we determined /^^ „, through numerical integration for each Uz, n^, n^, m, j 

and energy. We have found this approach to be computationally more time-consuming than the tabulation approach 
discussed above. 
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The evaluation of the integral , ^{£x',j), Eq. pO)) . proceeds analogously to that of , (ea'jJ)- We use [ioj 

T{a)U{aXx)^Y.TT^ (^8) 
— ' fc + a 

fc=o 

to separate out the energy-dependence that enters through Q^^ . The integral is then written as 

Iup,u'^,mi^x'J) = ^ lim^ ^ E T— HEH^'-^p'™' ^^^^ 

fc=0 2?7 

where 

Ci,n'^,n. - f Rnp,miP)RK,m (f) CXp (^-|^) L, (^^) pdp. (AlO) 

We tabulate C^^ ^ and calculate If^^ „, ,„(eA'jj) foi" ^ach e_^/ "on the fly". The convergence behavior of 
^np,n'p,mi^x'J) similar to that of I^.^n'^'^X' J)- 
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